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Abstract 



Inclusive momentum spectra are measured for all charged particles and 
for each of tt^, K^, K^/K^, and p/p in hadronic events produced via e"'"e~ 
annihilation at ^/s=58GeV . The measured spectra are compared with QCD 
predictions based on the modified leading log approximation(MLLA). The 
MLLA model reproduces the measured spectra well. The energy depen- 
dence of the peak positions of the spectra is studied by comparing the 
measurements with those at other energies. The energy dependence is also 
well described by the MLLA model. 



1 Introduction 

The LLA (leading log approximation) parton shower model well reproduces the 
various distributions of observables for hadronic final states of e~^e~ collisions, 
when the "coherence" effect of soft gluons is taken into account. Several Monte 
Carlo programs were written based on this scheme (JETSET63[]T]], for example) 
and were used in various experiments to study QCD. However, since the coherence 
effect is the consequence of higher order corrections, the effect could only be 
inserted by hand in the LLA scheme in these Monte Carlo programs. 

On the other hand, even in low region, the momentum spectrum of gluons 
in the parton shower process can be analytically calculated using the modified 
leading log approximation(MLLA) 0. The coherence effect is taken into account in 
MLLA by consistently importing a part of next-to-leading order corrections. The 
distribution of the particle spectra is expressed as a function of two parameters, 
F=log(^) and X=\ogm: 



r{B) 



cjA)^^(a,5 + i,cjA) (i; 



Here Xp is the momentum of a particle normalized by the beam energy E = ■\/s/2, 
9 is the opening angle of the jet cone, and Cp=^, 6 = y , and i? = |^ for 

five quark flavors, respectively. The two functions, $ and are two solutions of 
the confluent hyper-geometric equation. The quantities A and Qo are the QCD 
scale parameter and the energy cut-off of the parton evolution, respectively. 

This calculation predicts depletion of soft partons as a consequence of the 
destructive interference of soft gluons. This depletion shows up clearly in eq. |I| 
when written as a function of ^ = ln(l/xp). This function has a maximum at a 
certain ^ value and decreases in the larger ^ region. The coherence effect reduces 
the available phase space in this region and therefore the effect can be studied by 
comparing the measured inclusive cross section with this calculation. 

However, since this expression is calculated for partons, it cannot be compared 
directly with the measurement unless the distribution at the level of final state 
hadrons is ensured to be similar to that of partons. The concept of Local Parton 
Hadron Duality (LPHD) dictates that the distribution of the final state hadrons 
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is closely related to that of partons|P]. The conversion of partons into hadrons, 
which occurs at low virtuality scale, includes only small momentum transfer, and 
hence leaves the distribution essentially unchanged. In this article we compare the 
calculated parton spectrum (eq. |l]) directly with the measured hadron spectrum 
assuming LPHD. 

The Qo, which was primarily introduced to regularize coUinear singularity, 
gives a cut-off on parton energies. Therefore, the value of Qo should be close to 
the mass of the particle being considered [Q, while A should be common to all the 
particle species. These dependences can be experimentally tested by measuring 
the values of Qq and A for various particle species. 

The momentum spectrum might still be distorted by fragmentations and de- 
cays in spite of LPHD. This complication can be avoided partially when the energy 
evolution of the distribution is considered. The peak position of the distribution 
(eq. |1]) is given using the limiting spectra calculation in which A is assumed to be 
equal to QoH]: 



where is the number of colors (= 3). This should be compared with the 
variation as ^max ~ ^ expected from phase space consideration alone. 

2 Event Selection 

The data used in this analysis are accumulated by the TOPAZ detector at the 
TRISTAN e~^e~ collider at center-of-mass energies between 52.0 and 61.4 GeV. 
The average energy is 58.0 GeV, and the total integrated luminosity is 113.7/pb. 

Details of the TOPAZ detector is described elsewhere 0. In this analysis, the 
data from a Time Projection Chamber (TPC) is mainly used for tracking and 
dE/dx measurement for charged particles 0. The trigger condition relevant to 
this analysis is the track trigger which requires two or more charged tracks in the 
fiducial volume of the TPC, and the energy trigger which requires 4 GeV or more 
energy deposit in the barrel lead glass calorimeter. The event trigger is generated 
by a logical OR of those two, and the trigger efficiency for multihadronic events 
is practically 100%. 

Out of the triggered events, multihadronic events are selected by the follow- 
ing conditions; (a) five or more charged tracks having transverse momenta with 
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respect to the beam axis larger than 0.15 GeV/c are originating from the interac- 
tion point with the angle larger than 37° with respect to the beam axis, (b) total 
visible energy is larger than 1/2 of the center-of-mass energy, and (c) momentum 
imbalance along the beam direction is smaller than 0.4. By these conditions, the 
event selection efficiency is estimated to be 67.1% with a background contamina- 
tion of less than 2.0%. Two-photon processes and r+r^ productions are the main 
sources of the background. In addition to these cuts, the jet axis is required to 
have an angle larger than 40° with respect to the beam axis to ensure that the 
event is well contained in the detector acceptance. Applying these criteria, 11247 
events remain to be used in the analysis. 

3 Particle Identification and Measurement of Cross 
section 

The particle species vr^, and p/p are identified by measuring the dE/dx of 
tracks detected in TPC. The detail of the particle identification technique is de- 
scribed in ref. 1^. The typical resolution of the dE/dx measurement for minimum 
ionizing pions is 4.6%. Fig. |I| shows dE/dx distribution as a function of the track 
momentum for a part of the event sample. 

The cross section is calculated from the number of tracks in each momentum 
slice. In the low momentum region, the number of tracks for each particle species 
are directly obtained by counting the number of tracks in each dE/dx band. The 
dE/dx bands are, however, not well separated in the higher momentum region. To 
extract the number of each particle species in this region, the dE/dx distribution 
in a momentum slice is fitted with a superposition of four Gaussians corresponding 
to e^, TT^'s, i^^'s and p/p in each momentum slice. The widths and the centers 
of the Gaussians are determined from the measurement for Bhabha events and 
cosmic ray yu's. Only the normalizations of the Gaussians are the free parameters 
of the fit. A typical fit to the dE/dx distribution is shown in Fig. |^. Because of a 
hardware calibration problem of the TPC, the dE/dx resolution is time dependent. 
Therefore the data sample is divided into two groups and the procedure to count 
the number of tracks of each particle species is done separately for these two 
groups. 

Ks's are identified by searching for their daughter charged pion pairs in TPC. 
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Figure 1: dE/dx distribution as a function of track momentum measured by 
TOPAZ-TPC. 
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Figure 2: A typical fit to dE/dx distribution in a momentum bin. The distribution 
is fitted by four Gaussians corresponding to e"^, n"^, and p/p. 

A pair of tracks must satisfy following conditions to be identified as a Kg-, a) 
the distance of the closest approach of the two tracks is less than 0.8cm; b) the 
distance from the interaction point to the decay vertex is longer than 2.0cm where 
the decay vertex is defined as the center of the closest approach of two tracks; c) 
the angle formed by the vector from the interaction point to the decay vertex and 
the momentum sum of the two tracks at the decay vertex is less than 8 degrees if 
the position of decay vertex is longer than 6.0 cm; or c') the closest distance from 
the vector sum to the interaction point is less than 0.6 cm if the decay length is 
shorter than 6.0 cm; d) the closest approach of one of the tracks to the interaction 
point is more than 0.5cm if the momentum sum of the tracks is less than 2GeV/c; 
e) tracks are identified as vr^ by the dE/dx measurement in the TPC; and f) the 
tracks are not from the gamma conversion. 

For each of track pairs remained after these selections, the invariant mass is 
calculated assuming the tracks are charged pions. Fig. |] shows the mass distribu- 
tion for the reconstructed Ks's. The distribution is fitted by a sum of a Gaussian 
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M(7T^7T ) (GeV/c^) 



Figure 3: Reconstructed mass distributions of Kg in two different momentum bins: 
(a) 1.8< i <2.1, (b) 3.3< i <3.6 

and a background function (exponential + polynomial) and the number of recon- 
structed KsS is obtained to be 771 ± 61. In the same way, the numbers of these 
KgS are counted for each momentum slice. 

The inclusive cross section, l/(Jhad dcr/dxp is then calculated from the number 
of particles observed in each momentum slice using following formula: 

-^^-A{x.)^^^^, (3) 

ahad dXi AXi Nhad 

where Xi is the i'th slice of the momentum fraction, ahad is the total hadronic 
cross section, Axi is the width of the slice, N^ad is the number of hadronic events 
in the event sample, and Nobs{xi) is the number of particles counted in the slice, 
respectively. A{xi) is the factor to correct the counted number for the effects 
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of detector acceptance and initial state radiation. A{xi) is calculated for each 
momentum slice separately with the Monte Carlo simulation as 

^y-^'^J ATtotal I KTtotal 

gen sim 

where Ngen{xi) is the generated number of particles in the slice Xi without includ- 
ing the effects of the detector acceptance and the initial state radiation. N^f^'' is 
the corresponding total number of generated events. Nsim{xi) is the number of 
detected particles in the slice Xi generated with the effects of detector acceptance 
and initial state radiation, and A^*°^' is the corresponding total number of gen- 
erated events. The JETSET6.3|]T]] and 7.3||^ Monte Carlo programs are used to 
obtain these numbers combined with the TOPAZ detector simulation program. 
As for Kg, to convert the counted number to the cross section of K^/K'^, A{xi) is 
multiplied by 2. 

The cross sections are measured as a function of C,- The measured cross section 
for all charged particles is shown in Table |l|. Table shows the cross sections 
measured for each of vr^, K^, p/p, while that for / K'^ is shown in Table 0. 

The systematic error in the measurement is estimated by considering following 
sources. One is the systematic ambiguity in the acceptance correction. This is 
studied by varying the fragmentation parameters in the Monte Carlo programs 
and by turning on and off the energy loss and the nuclear interaction effects in the 
detector simulation program. The ambiguity is estimated to be 2% for all charged 
particles while 3% for vr^, K"^ and p/p. For vr^, and p/p, the ambiguity in the 
dE/dx curve as a function of momentum used to extract numbers of each particle 
species is also a source of the systematic error. This is estimated to be 5% for vr^ 
and 10% for and p/p. 



4 Comparison with MLLA prediction 

The measured particle spectra are then compared directly with the MLLA pre- 
dictions by assuming LPHD. Since the A and Qq in the MLLA calculation are 
unknown, their values are determined by fitting the MLLA formula of the cross 
section to the data. Fig. § shows the measured cross section for all charged parti- 
cles as a function of ^ together with the result of the best fit shown in the solid line. 
As seen from Fig. |, the fitted MLLA calculation well reproduces the measured 
cross section. 
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Table 1: Measured cross section for all charged particles as a function of ^. Errors 
include both of statistical and systematic errors. 
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Table 2: Measured cross sections for tt^, and p/p as functions of ^. Special 
non-uniform binning of ^ is used to optimize the particle identification by the 
dE/dx measurement. Errors include both of statistical and systematic errors. 
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^ = ln(l/xp) 


l/cTtoi dcr{Ks)/d^ 


1.80 


0.58±0.10 


2.25 


0.65±0.12 


2.55 


0.68±0.10 


2.85 


0.61±0.09 


3.15 


0.62±0.09 


3.45 


0.49±0.09 


3.80 




4.25 




4.05 


0.23±0.04 



Table 3: Measured cross section for K^/ as a function of C,- Errors are statistical 
only. 

Fig. ^ shows the cross sections measured for vr^, K^, K^/K^ and p/p with the 
fitted MLLA calculations. Also shown are the measurements by PEP4/TPC0 
at a/s = 29GeV for vr^, K^, and p/p. These measurements are also fitted by 
the MLLA calculations. The fit is performed in the range where the numerical 
calculation of the MLLA function is reliable. The range is typically 1.0< C, <4.0. 
The measured cross sections are well reproduced by MLLA. The figure demon- 
strates the energy evolution of the peak position of the spectra: the peak position 
in ^ becomes larger in TOPAZ measurements at ^/s = 58GeV than those of 
PEP4/TPC at y/s = 29GeV regardless of particle species. This is also predicted 
by the MLLA calculation. The detail of the study of the evolution of the peak 
position is described in the next section. 

Tables ^ and ^ summarize the obtained values of A and Qq. Two different 
fits are performed: one is to determine both of A and Qo from the fit, and the 
other is to determine Qq from the fit with A fixed at 200 MeV. In both cases, the 
value of Qo becomes larger for heavier particles. Meanwhile, in the first case, the 
determined value of A stays at relatively lower values (100-300 MeV). This result 
gives a strong support to the MLLA + LPHD conjecture. 
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TOPAZ(58GeV) 


PEP4(29GeV) 


Free fit 


A (MeV) 


(5o(MeV) 


A (MeV) 


go(MeV) 


all charged 


291±10 


375±8 






7r± 


281±20 


339±25 


270±7 


250±7 


X± 


118±68 


575±80 


243±74 


531±80 




185±19 


649±76 






viv 


143±81 


657±90 


356±105 


531±100 



Tabic 4: A and Qo for each of particle species determined from the fits with both 
A and Qo as free parameters. 





TOPAZ(58GeV) 


PEP4(29GeV) 


A = 200McV 


Qo (McV) 


Qo (McV) 


all charged 


297±12 






275±12 


217±6 


X± 


588±44 


514±45 


K^W 


663±110 




P/P 


627±92 


390± 50 



Tabic 5: Qq for each of particle species determined by the fits with A fixed at 
200MeV. 
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Figure 4: The cross section measured for all charged particles as a function of 
^ = ln{l/xp). Solid curve shows the fitted MLLA calculation. 

5 The Energy Evolution of ^ Distribution 

The energy evolution of the peak positions in ^ is studied by comparing the mea- 
surements at yi=29GeV(PEP4/TPC), 58GeV(TOPAZ) and 91GeV(ALEPH)|0| 
as a part of the PTA project Table ^ shows the result. The peak positions 
and their errors are estimated using the fitted MLLA functions for PEP4 and 
TOPAZ measurements. The MLLA predictions using the limiting spectra calcu- 
lation (eq. ^ with A = 200MeV are also shown. As seen, the peak positions 
for light particles (all charged particles and vr^) agree with the limiting spectra 
calculation. The peak positions become smaller for heavier particles regardless of 
the energy scale. 

From eq. ^ the MLLA prediction of peak position can be approximated in the 
linear form alni^^ + b where a to be 0.5. If the gluon coherence effect is not taken 
into account (phase space only), a should become 1.0. To obtain the value of a 
from the measured peak positions, a linear fit is performed to the peak positions 



12 




Figure 5: The cross sections measured for (a) n^, (b) K^, (c) p/p and (d) /K^. 
Both of TOPAZ and PEP4/TPC measurements are plotted with MLLA fits (Sohd 
curve:TOPAZ, Dashed curve:PEP4). 
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PEP4/TPC(29GeV) 


TOPAZ(58GeV) 


ALEPH(91GeV) 


all charged 




3.30±0.05 


3.61±0.01 




3.10±0.12 


3.48± 0.18 


3.81±0.02 




2.30±0.20 


2.56± 0.29 


2.63±0.04 






2.90±0.24 


2.88±0.03 


viv 


2.30±0.20 


2.50±0.29 


3.00±0.09 


MLLA(A = 200MeV) 


3.02 


3.46 


3.74 



Table 6: The energy evolution of the peak positions in ^. Also shown are predic- 
tions obtained using the limiting spectra calculation of MLLA. 



as a function of a/s for all charged particles measured by TASS0[|12|, TOPAZ, 
and ALEPH as shown in Fig. ^ From the fit, the value of a is obtained to be 
0.54 ± 0.05. A similar fit is also done to the peak positions for charged pions 
(TASSO, PEP4/TPC, TOPAZ and ALEPH) and a is measured to be 0.58 ±0.04. 
These values are consistent with the MLLA predictions. 



6 Summary 

The inclusive cross sections are measured as a function of ^ = ln{l/xp) for all 
charged particles and for each of vr^, K^, p/p and / K'^ in the hadronic events 
taken at y^=58GeV. The cross sections were compared with the MLLA calcu- 
lation assuming LPHD. The MLLA formula describes the observed distributions 
very well for all particle species over the wide momentum range. 

By this comparison, Qq and A in the MLLA expression are determined for 
each particle species. The determined Qq for each of the particle species is close 
to the mass of the particle, while the A is almost constant for these particles. 
The Qo values are also determined with A fixed at 200MeV. The obtained Qq also 
coincides with the mass of each particle. The same analysis is carried out for the 
measurements at ^/s = 29GeV by PEP4/TPC and similar tendency is observed. 
These results show that the momentum distribution of particles is identical to 
that of partons at the end-point of the parton shower evolution independently of 
the energy scale. This supports the MLLA + LPHD conjecture. 

The energy evolution of the peak position in da/d^ is studied by comparing 
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Figure 6: The peak position in ^ as a function of the square of the center-of-mass 
energy. Sohd hne shows the hnear fit to the peak positions of charged particles 
while dashed line to those of charged pions. 

our data with the measurements by PEP4/TPC and ALEPH. The measured peak 
positions of light particles are well reproduced by the limiting spectra calculation. 
The measured peak positions arc fitted to a linear function of log of the beam 
energy with the measurements at other energies. The obtained slope of the linear 
function is consistent with the MLLA prediction, while it excludes models without 
the gluon coherence effect. 
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